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The  strictly  mathematical  ocean  tide  model  developed  in  Part  I 
of  this  paper  is  modified  in  order  to  include  realistic  hydrody¬ 
namical  barrier  effects  of  narrow  ocean  ridges  and  other  large 
bottom  irregularities.  This  modification  begins  with  a  hydro- 
dynamical  redefinition  of  the  ocean  bathymetry  at  over  3,000 
grid  points,  increasing  simultaneously  the  depth  data  range  to: 
10m  -*  7,000m.  In  a  second  step  a  unique  hydrodynamical  in¬ 
terpolation  technique  is  developed  that  incorporates  into  the 
model  over  2,000  empirical  tide  data  collected  around  the 
world  at  continental  and  island  stations.  This  interpolation  is 
accomplished  by  a  controlled  cell-wise  adjustment  of  the  bottom 
friction  coefficient  and  by  allowing  a  monitored  in-  or  out-flow 
across  the  mathematical  ocean  boundary  and  so,  redefining  a 
more  physical  shoreline.  Extensive  computer  experiments  were 
conducted  to  study  the  characteristics  of  the  novel  friction  laws 
and  hydrodynamical  interpolation  methods.  The  computed  Mj- 
tide  data  along  with  all  (specially  labeled)  empirical  constants 
are  tabulated  in  map  form  for  four  typical  30°  by  50°  ocean 
areas.  It  is  estimated  that  the  tabulated  tidal  charts  permit. a 
prediction  of  the  Mt-tide  elevation  of  the  ocean  surface  over 
the  geoidal  level  with  an  accuracy  of  better  than  5  cm  anywhere 
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in  the  open  ocean  and  with  somewhat  less  accuracy  near  rough 
shorelines.  With  the  forthcoming  construction  of  the  lesser  S2, 
N2,  and  K2;  Kt,  Ot,  P,,  and  Qi;  and  Mf,  Mm,  and  Ssa  tidal  con¬ 
stituents,  the  total  tide-prediction  error  can  be  kept  below  the 
10-cm  bound  posed  by  applied  researchers  of  today. 

In  Part  I  of  this  paper  (Schwiderski,  1979)  a  purely  hydrody- 
namical  ocean  tide  model  has  been  developed  and  tested.  This 
model  has  been  applied  to  compute  a  preliminary  M2  ocean  tidal 
chart  (Schwiderski,  1976, 1).  (References  listed  in  Part  I  are  indi¬ 
cated  by  the  added,  I,  after  the  year  specification.)  The  results 
were  found  encouraging  and  satisfactory  for  some  applications. 
However,  significant  shortcomings  still  persisted,  especially  over 
narrow  ocean  ridges.  The  remaining  deficiencies  were  attributed  to 
local  distortions  and  retardations  of  the  tidal  waves  due  to  hydro- 
dynamical  barrier  effects  of  ocean  ridges  and  other  bottom  ir¬ 
regularities. 

In  the  following  Part  II  of  the  paper,  an  attempt  will  be  made  to 
eliminate  the  shortcomings  of  the  purely  theoretical  model  by  using 
a  hydrodynamically  defined  bathymetry  of  the  oceans  and  by  in¬ 
corporating  directly  empirically  known  tide  data  into  the  discrete 
tide  model  described  in  Part  I.  The  latter  modification  will  be  ac¬ 
complished  by  a  controlled  local  adjustment  of  the  bottom  friction 
coefficient  and  by  allowing  a  monitored  in-  or  out-flow  across  the 
mathematical  ocean  boundary,  and  thus  redefining  implicitly  a 
more  physical  shoreline.  A  detailed  discussion  of  the  quality  of  the 
new  Ms  ocean  tidal  charts  will  be  given  in  the  sections,  “Quality  of 
the  Ocean-Tide  Model”  and  “Conclusions.” 

The  complete  Ms  ocean  tide  is  published  in  tabulated  map  form 
in  Schwidenki  (1979c  I).  Similar  charts  for  the  Sa,  Ns,  Ks,  Oi,  Pi, 
Qt,  Mf,  Mm,  and  Ssa  ocean  tides  (see  Part  I,  Table  1)  are  under 
construction  and  will  be  published  in  additional  papers.  A  separate 
tabulation  of  the  new  hydrodynamically  defined  ocean  bathymetry 
may  be  found  in  Schwiderski  ( 1978a  I).  All  tidal  and  bathymetry 
data  will  be  available  in  tape  form  at  the  Naval  Surface  Weapons 
Center,  Dahlgren,  Virginia  22448. 
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Hydrodynamically  Defined  Ocean  Bathymetry 
A  reinspection  of  the  bathymetric  data  revealed  clearly  that  even  a 
1*  by  1°  grid  scheme  falls  far  short  in  representing  a  narrow  ocean 
ridge  in  a  hydrodynamically  proper  fashion.  The  defect  is  partic¬ 
ularly  compounded  when  the  narrow  ridge  parallels  a  deep  trench. 
The  reason  for  this  deficiency  is  obviously  the  purely  hydrostatic 
character  of  the  averaging  principles  employed  by  Smith  et  al. 
(1966  I)  in  order  to  assign  a  depth  value  at  the  center  of  a  mesh 
cell  that  is  supposed  to  be  representative  for  the  entire  cell.  For 
instance,  if  the  area  of  a  mesh  cell  is  (by  subjective  sight)  more 
than  half  land,  then  it  is  called  a  “land  cell,"  and  the  cell  is  given 
(for  the  present  purpose)  the  depth  value  “zero.”  In  the  alterna¬ 
tive  case,  the  cell  is  declared  “oceanic,"  and  a  depth  value  is 
assigned  that  conserves  the  estimated  actual  water  mass.  Because 
of  those  hydrostatic  principles,  cells  were  found  that  contained 
elongated  islands  crossing  even  several  cells,  but  every  cell  was 
declared  oceanic.  Moreover,  an  oceanic  trench  portion  of  die  cell 
with  some  7,000-m  true  depth  produced  an  average  depth  of  more 
than  3,500  m.  Clearly,  for  ocean  current  models  the  entire  cell  rep¬ 
resents  an  impassable  wall,  and  the  depth  value  should  be  “zero" 
instead  of  3,500  m. 

hi  order  to  eliminate  the  shortcomings  of  the  bathymetric  data 
compiled  by  Smith  et  al.  (1966  I),  the  depth  values  were  revised 
by  using  the  following  “hydrodynamical"  principles: 

(a)  Boundary  cells  at  or  near  continental  shorelines  con¬ 
sisting  of  more  than  half  oceanic  areas  of  depths  larger  than 
5  m  were  designated  ocean  cells,  and  the  average  oceanic 
depth  values  were  assigned  as  die  “hydrodynamically"  aver¬ 
aged!  depths  to  the  entire  cells.  The  new  depth  value  is  pref¬ 
erable  to  the  "hydrostatically"  averaged  depth,  which  pre¬ 
serves  the  actual  water  mass  but  ascribes  artificially  a  shallow 
shelf  character  to  the  cell. 

(b)  Maud  calls  were  declared  terrestrial  oeBs  with  depths 
aero  if  either  the  island  areas  were  larger  than  half  die  mesh 
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areas  or  the  (elongated)  island  lengths  exceeded  the  mesh 
diameters. 

(c)  Island  cells  that  remained  oceanic  cells  were  assigned 
depth  values  less  than  the  hydrostatically  averaged  values.  In 
this  case  and  in  situations  of  submerged  seamounts  or  narrow 
ocean  ridges  (e.g.,  Aleutian,  Marianas,  and  Caribbean),  the 
hydrodynamical  depths  depended  on  the  assessed  “barrier’' 
effects  of  the  obstacles:  the  longer  and/or  higher  the  barrier, 
the  lesser  the  depth.  In  general,  the  average  “ridge  depth’’  was 
assigned  to  the  entire  cell. 

(d)  The  assigned  minimum  depth  (Part  I,  Equation  50)  was 
lowered  to 

Hm  =  min  H(A,0)  =  20  m,  (1) 

which  is  further  lowered  to  10  m  by  the  averaging  Equations 
65  in  Part  I.  (All  notations  of  Part  I  are  used  unchanged  in 
the  present  paper.) 

The  hydrodynamically  justified  principles  (a)  to  (c)  are,  natu¬ 
rally,  quite  subjective  and  by  no  means  free  of  any  error.  Never¬ 
theless,  some  computational  experiments  indicated  only  very  minor 
effects  of  isolated  depth  data  changes.  More  than  3,000  depth 
values  were  changed,  but  only  very  few  of  those  required  addi¬ 
tional  readjustments  in  order  to  keep  some  limitation  on  the  first 
and  second  derivatives  of  H(x,6);  i.e.,  on  the  relative  differences 
given  by  Equation  65  in  Part  I.  Furthermore,  the  hydrodynamical 
interpolation  of  empirical  tidal  data  (section,  “Hydrodynamical 
Interpolation  of  Empirical  Tide  Data’’)  known  at  continental  and 
island  stations  greatly  diminishes  the  need  for  precise  boundary- 
depth  data.  The  revised  depth  data  bank  used  in  the  new  tidal 
computations  are  published  in  Schwidenki  (1978a  I). 

Empirical  Tide  Data 

The  new  tide  model  incorporates,  by  a  unique  hydrodynamical  in¬ 
terpolation  procedure  (next  section),  empirical  tidal  data  observed 
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and  harmonically  analyzed  at  numerous  continental  and  island 
stations.  These  data  were  taken  from  publications  by  the  National 
Ocean  Survey  (1942),  the  International  Hydrographic  Bureau 
(1966),  the  British  Admiralty  ( 1977 1),  and  by  Pekeris  and  Accad 
(19691),  Zahel  (19701, 1973 1),  Cartwright  (1971),  and  Luther 
and  Wunsch  (1975  I).  Unfortunately,  the  most  recent  publication 
by  the  British  Admiralty  lists  harmonic  constants  only  for  die  four 
major  tide  components  M2,  S2,  Ki,  and  Oi  and  excludes  the  Euro¬ 
pean  waters  completely. 

'  The  voluminous  data  banks  had  to  be  screened  in  order  to  elimi¬ 
nate  observations  that  are  meaningless  or  unreliable  for  the  present 
ocean-tide  investigations.  For  example,  tidal  constants  were  ex¬ 
cluded  that  were  listed  for  stations  deep  inside  estuaries  or  harrow 
bays  (eg.,  Hudson  River,  Bay  of  Fundy),  at  the  mouths  of  large 
rivers  (e.g.,  Amazon),  between  sheltering  islands  (e.g.,  Alexander 
Archipelago,  Solomon  Islands),  and  inside  sheltering  reefs  (e.g., 
Great  Barrier  Reef) . 

About  2,500  stations  were  selected  for  further  examination  of  their 
data  concerning  locally  restricted  distortions.  For  instance,  some  data 
taken  over  short  distances  along  a  coastline  displayed  rather  drasti¬ 
cally  alternating  times  of  high  water,  which  are  obviously  meaning¬ 
less  for  oceanic  tidal  studies.  At  many  stations,  different  tables  give 
different  tidal  constants.  Some  of  those  discrepancies  at  island  sta¬ 
tions  are  shown  for  the  Ms-tide  in  Table  1.  Similarly,  for  some  mesh 
cells,  several  different  station  data  were  available,  and  only  one 
representative  average  had  to  be  chosen.  This  situation  is  illustrated 
in  Table  2  for  the  Mr-tide  around  Bermuda..  Many  of  those  differ¬ 
ences  can  probably  be  explained  as  simple  errors  in  printing  or 
computing.  For  instance,  the  phase  difference  of  about  1  hr  at 
Port  Galets  on  La  Reunion  Island  (Table  1)  seems  to  be  due  to 
some  error  in  observing  die  correct  reference  time,  which  varies 
from  listing  to  listing.  Most  differences,  such  as  those  shown  for 
^bnush  stations  in  Table  2,  are  definitely  true  local  variations.  In 
'fMa  connection,  the  important  tidal  measurements  by  Gafiagher  et 
d.  (1971)  at  Fanning  Atoll  in  the  central  Pacific  any  be  men- 
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Tabic  1 

_ Empirical  M,-tkta  differences. _ 

Station  B.A.T.(77)>  N.O.SJ42, ?  Others  Initialed 


Latitude,  Longitude 

P(m)  tH°) 

Hm) 

«r°; 

SfO) 

Tenerife,  Canary  island 

(A) 

0.87 

0.89 

Z* 

28*2974  16°14'W 

18 

30 

Port  Praia,  Cape  Verde  i. 

(A) 

0.42 

0.43 

Z* 

14#55'N  23°31'W 

244 

220 

Ascension  Island 

(A) 

0.33 

0.51 

p,*z 

7#55'S  14*25^ 

177 

174 

St.  Helena  Island 

(A) 

0.32 

0.34 

PX 

15°55'S  5°42'W 

81 

87 

Tristan  da  Cunba  Island 

(A) 

0.23 

0.34 

px 

37#02/S  12°18'W 

12 

354 

Agalega  Island 

(1) 

0.29 

0.29 

z 

10*26*8  56°40*E 

350 

290 

Port  Qalets,  La  Reunion  1. 

(1) 

0.16 

0.14 

0.14 

z 

20*55'S  56°17'E 

302 

328 

328 

Mawson,  Antarctica 

(1) 

0.04 

0.04 

z 

87°38'S  82°53'E 

232 

155 

Wilkes  Station,  Antarctica 

(1) 

0.28 

0.36 

z 

avis's  110°31*E 

162 

140 

Welles  Harbor,  Midway  1. 

<P) 

0.11 

0.11 

PX 

28*12*N  117*22*W 

62 

91 

Eniwetok  Atoll,  Marshall  l. 

<P> 

0.36 

0.36 

11*21*N  162°21'E 

127 

137 

Nee  Wallis.  Fiji  Island 

(P) 

0.53 

0.52 

13*22*8  178*11*W 

178 

154 

Suva  Harbor,  Vrti  Levu, 

FIJI  (stand 

(P) 

0.56 

0.50 

18° 00*8  178°28'E 

196 

212 

•  B.A.Tj77)j«BtmNia  AdfflfrUy  TaMM  (1»77  I). 

•»  »  Udal  pSaaa  fitttv  to  Owmrtca. 
anas.  HO  «  Naaowl  Oaaan  kmr 
*Z>2M  (1t70  I). 

/Pat  Catorte  and  Aoead  (INC  I). 

dosed.  Tides  outride  and  inside  the  small  atoll’s  Ippm  dMfcnd  by 
•boot  30%  (20  on)  in  — lptitoda  mad  by  a  phase  lag  of  iboot  30* 
(1  hr,  40  ado.). 

In  general,  the  moct  recent  listings  in  the  British  Admiralty  Tide 
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Tildes  (1977  I)  were  chosen  over  older  tabulations  as  the  most 
reliable.  The  selection  of  the  data  was  further  aided  by  earlier  and 
subsequent  tidal  computations.  Altogether,  some  1,700  Mj-tide 
data  were  selected  and  assigned  to  the  craters  of  their  respective 
mesh  cells.  Using  linear  interpolation  and  tidal  computations,  the 
total  number  of  prescribed  tide  data  used  in  the  Ms-tide  construc¬ 
tion  was  increased  to  more  than  2,000.  Essentially  all  continental 
boundary  cells  carry  empirically  supported  tide  data.  The  empirical 
coverage  is  only  marginal  at  arctic  and  antarctic  shorelines.  Most 
empirical  tide  data  known  at  island  stations  are  also  included  in  the 
tide  model. 

Naturally,  it  must  be  remembered  that  the  selection  of  represen¬ 
tative,  empirical  tidal  data  (compare  depth  data,  section  before 
MHydrodynamically  Defined  Ocean  Bathymetry*)  is  not  at  all  free 
of  subjective  judgment  and  may  be  somewhat  erratic.  Obviously, 
only  future  additional  tidal  measurements  can  improve  this  model 


Station 

Latitudo,  Longitude 

(•{cm) 

nfffffrrCv 

8t  George’s  Island 
a2JSN.04.7BW 

36 

360 

British  Admiralty  (19771) 

8t  David’s  Island 
32J7N.04J6W 

34 

366 

British  Admiralty  (10771) 

Orest  Sound 
32J2N.04JSW 

at 

6 

British  Admiralty  (1977 1) 

8L  George's  Island 
aSJ7N.04.7OW 

as 

0 

National  Ocean  Survey  (1942) 

of.  utoryi  t  HM9KI 
08.4ON.04.7OW 

37 

0 

Pekoris  and  Ac cad  (1909 1) 

•i  George's  Island 
0SJ7N.04.7OW 

as 

360 

Zahsl  (1970 1) 

•t  George's  Island 
aS.4ON.04.7OW 

as 

366 

Zettler  at  al.  (1979) 

Deep  Sea  (GOBI  IV) 
SSJ0N.O4.SOW 

aa 

11 

J.  T.  Kuo  Letter  (1977) 
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in  this  respect  Nevertheless,  according  to  the  instruction  notes  ac¬ 
companying  the  British  Admiralty  Tide  Tables  (1977  I),  it  can 
probably  be  assumed  that  almost  all  important  tide  data  selected 
cany  an  accuracy  that  is  at  least  as  high  as  the  desired  10  cm.  In 
any  case,  computational  experiments  showed  that  isolated  reason¬ 
able  variations  of  the  boundary-tide  data  do  not  affect  significantly 
the  adjacent  oceanic  tides.  It  was  also  found  insignificant  to  the 
overall  quality  of  the  tide  model  whether  the  empirical  data  were 
assigned  to  the  centers  or  to  the  shore  boundaries  of  die  respective 
cells. 

Attempts  were  made  to  incorporate  also  recent  deep-sea  tidal 
measurements  into  the  present  model.  Since  die  hydrodynamical 
interpolation  of  empirical  data  is  essentially  based  on  bottom  and 
boundary  irregularities  (see  next  section  (a)-(d)),  no  physically 
valid  justification  was  found  to  include  distant  offshore  deep-sea 
measurements  into  the  model.  However,  some  deep-sea  measure¬ 
ments  near  rough  shore  and  bottom  areas  were  included.  Fortu¬ 
nately,  without  exception,  all  excluded  offshore  deep-sea  measure¬ 
ments  known  to  the  author  agree  very  well  with  the  computed  M2- 
tide  data  (see  Table  3 ) . 


Table  3a 

Peap  aaa  Ma-tMe  data  lor  the  Quit  of  Mexico  and  Caribbean  Sea. 


Station 

Observed 

Model 

Error 

Latitude,  Longitude 

P(cm) 

H cm ) 

MV 

Afifcm) 

*MV 

W.  Florida  Shelf  St. 

7 

7 

0 

26.71  N,  84.25W 

Deep  Quit  St. 

1.3 

97 

1.6 

92 

+0.3 

-5 

24.77N,  89.65W 

Mister)  osa  Bank 

8 

226 

9 

22S 

+1 

-1 

18.88N,  83.81  W 

Rosalind  Bank 

7 

84 

8 

89 

+1 

+5 

16.61  N,  80.34W 

East  Carib.  St  (6-month) 

0.5 

107 

1.6 

102 

+1 

-5 

16.54N,  64.6SW 

East  Carib.  St  (1-month) 

0.6 

156 

1.5 

151 

0.9 

-5 

16.52N,  64.91  W 

153 

146 

-5 

•  4  =  tkM  amplitude. 

»  «  a  tMal  phiaa  ralatlva  to  Qraanwleti. 
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Table  3b 

Deep-sea  M,-tide  data  for  the  Pacific  and  Atlantic  oceans. 


Station 

Observed 

Model 

Error 

Latitude,  Longitude 

i°(cm)  sb(°) 

Horn) 

H°) 

Atfcm)  m(°) 

Pacific  St.  1  (Middleton) 
58.76N,  145.71W 

110 

284 

included 

Pacific  St.  3  (Tofino) 
48.97N,  127.29W 

99 

239 

included 

Pacific  St.  (San  Francisco) 
38.16N,  1 24.91  W 

54 

227 

included 

Pacific  St.  (Josie  II) 
34.00N,  144.99W 

27 

267 

27 

273 

0 

+6 

Pacific  St.  (Fticki) 

32.24N,  120.86W 

43 

149 

included 

Pacific  St.  (Josie  1) 

31.03N,  119.80W 

43 

142 

included 

Pacific  St.  (Kathy) 

27.75N,  124.37W 

29 

128 

27 

130 

-2 

+2 

Pacific  St.  (Filloux) 

24.78N,  129.02W 

19 

107 

18 

105 

-1 

—2 

Atlantic  St.  1  (N.Y.  Bight) 
39.32N,  64.36W 

44 

350 

included 

Atlantic  St.  (N.C.  St.  1) 

48 

46 

—2 

32.69N,  75.62W 

Atlantic  St  (Savannah  B) 

88 

356 

included 

358 

+2 

31 .9514,  80.68W 

Atlantic  St.  (Scope) 

45 

15 

46 

+2 

30.43N,  76.42W 

Atlantic  St.  (AOML  1) 

34 

358 

35 

3 

+1 

+5 

28.1 4N,  69.75 W 

Atlantic  St  (AOML  3) 

34 

1 

34 

6 

0 

+5 

28.24N,  67.54W 

Atlantic  St.  (MERT) 

34 

359 

34 

4 

0 

+5 

27.99N,  69.67W 

Atlantic  St.  (REIKO) 

35 

360 

34 

6 

-1 

+6 

27.97N,  89.87W 

Atlantic  St  (EOIE-May) 

32 

1 

32 

6 

0 

+5 

28.48N,  89.33W 

Atlantic  St.  (EDIE-March) 

31 

3 

32 

7 

+1 

+4 

26.45N,  69.32W 

1 

7 

+6 

•  |  =  tidal  amplitude. 

*»  =  tidal  phase  relative  to  Greenwich. 
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Of  course,  the  continuity  gap  (Equation  4)  can  be  attributed  to 
the  following  major  causes  which  are  physically  plausible: 

(a)  The  bottom-friction  coefficient,  b  (in  A4  and  B4  of  Equa¬ 
tions  62  in  Part  I),  which  is  most  effective  in  boundary  cells, 
depends  on  local  shore  features  such  as  true  cell  size  and  bot¬ 
tom  slope  and  roughness. 

(b)  The  boundary  cells  are  idealized  by  definition  of  strictly 
mathematical  boundaries  (see  Figure  1 ). 

(e)  The  depth  data  of  boundary  cells  are  subjectively  de¬ 
fined  and,  hence,  faulty  (section,  “Hydrodynamically  Defined 
Ocean  Bathymetry” ) . 

(d)  The  empirical  tidal  constants  in  Equation  3  are  also 
faulty  to  some  degree  because  of  inaccurate  measurements, 
harmonic  analyses,  and  subjective  selections  and  assignments 
to  the  centers  of  the  boundary  cells  (preceding  section). 

(e)  The  discrete  ocean-tide  model  is  certainly  not  an  exact 
description  of  the  true  oceanic  tide;  e.g.,  at  boundaries,  non¬ 
linear  inertial  terms  assume  significance. 

Obviously,  the  last  two  (hopefully  minor)  faults  can  be  reduced 
only  through  continued  future  observations  and  modeling.  How¬ 
ever,  the  first  two  faults,  (a)  and  (b),  can  be  weakened  by  “hydro¬ 
dynamically  interpolating”  the  empirical  tidal  elevations  (Equation 
3)  into  the  tidal  model  and  narrowing  the  continuity  gap  (Equa¬ 
tion  4)  to  an  acceptable  level  as  follows: 

(A)  Adjusting  the  velocity  field  by  a  locally  controlled  im¬ 
plicit  variation  of  the  bottom-friction  coefficient,  b,  in  Equa¬ 
tions  62  Part  I. 

(B)  lifting  the  strict  condition  of  no-flow  across  the  mathe¬ 
matical  ocean  boundary  and  allowing  for  a  monitored  in-  or 
out-flow  by  implicitly  defining  a  more  physical  ocean  boun¬ 
dary  (Figure  1). 

As  was  pointed  out  in  Part  I,  section.  The  Discrete  Ocean-Tide 
Equations  (DOTE*),”  due  to  the  choice  of  the  finite-difference 
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TRUE  BOUNDARY 
MATH.  BOUNDARY 


(b) 


figure  1.  Boundary  cell  in-  and  out-flow  illustration:  (a)  and  (i),  also  (b) 
and  (E)  are  half -periods  apart.  (Shaded  region  is  land  area.) 

parameter  *  =  1,  the  bottom  friction  coefficient,  b,  in  A4  and  B4  of 
the  momentum  equations  (Part  I,  Equations  60)  can  be  considered 
implicitly  varied  in  the  mesh  cell  &•,>  by  directly  replacing  the 
velocity  components  in  Equation  2  as  follows: 
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and 


v!£  -» +  |i^j|  (wm + wbi), 
"S,,.  -  £C,‘,,  +  IdJ  <w«.  +  wb,), 
VIS  -*  V£i  +  |f»«|  (wVl  +  WV,), 

d,  -  d, + id,  i  <wv« +«*•), 


(5) 


provided  ¥•  0;  i.e.t  provided  an  empirical  tidal  amplitude  is 
available  for  the  considered  mesh  cell.  In  Equation  5,  the  consis¬ 
tency  and  scale  parameters  (u,U)  and  (v,V)  are  defined  by 


mid 


Ui  =  1,  Di  =  0 

forA{S,l-tfff‘<0, 

Mi  =  0,  Ql  =  ilii* 

otherwise,  but 

Mr 

(6a) 

Ml  =  0,  Bl  =  0 

ifS-^O; 

Ml  =  1,  Hi  =  0 

for  •  !/&?*»  >  0, 

(6b) 

M«  88  0,  Bs  88  Am+pn 

otherwise; 

© 

II 

£ 

* 

II 

£ 

fOTA^;.^i<0, 

Vi  =  0,  Vl  =  Bit,* 

(6c) 

otherwise; 

* 

II 

* 

II 

© 

for  >  0, 

Vs  =  0,  Vi  =  Bl,*-! 

otherwise,  but 

(6d) 

Vs 38  0,  Vs  88  0 

The  continuity  gap  (Equation  4)  will  be  narrowed  when  the  “con¬ 
trol  parameters”  w  and  W  are  determined  successively  by: 


w1 


10 

with  tiie  tint  “control  limit” 


for  C  *  0, 
for  {  =  0 


(7a) 


1 
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and 


w  j  <  k\ 


[A^i-  wt]/Z 
0 


for  i  *  0, 


for  £  =  0 


with  the  second  control  limit 


w\  ^  ki 


(7b) 

(8a) 


(8b) 


where  (see  Equation  3) 


t = Cl  [u,  IUBI  +  Im«  +  v,ci|K»;i  +  V,C*J  V£-J_,h 
and  M9) 

ci  [si  iicii  rasj]  +v2c«|f^_j) 

It  is  important  to  note  that  ui  •Hi  =  Oand  Vi  •  "Vi  =  0  for  i  =  1,2. 
Accordingly,  both  control  limits,  ki  and  k2,  which  are  at  one’s  dis¬ 
posal,  regulate  the  allowed  decrease  or,  respectively,  increase  of  the 
velocity  components  in  Equations  5 ;  i.e.,  the  implicitly  permitted  cor¬ 
responding  increase  or  decrease  of  the  local  bottom-friction  coeffi¬ 
cients.  Since  the  integration  sweeps  across  the  ocean  from  m  “  p  to 
360  and  n  =  4  to  168,  the  special  choice  of  ui  =  Hi  =  0  and  v2  ='V2 
=  0  in  Equations  6a  and  6d  excludes  possible  double  adjustments 
of  the  velocity  components.  Also,  if  ut  Th  and/or  v2  **  tfe,  back¬ 
ward  adjustments  of  the  tidal  elevations  via  the  corresponding 
Equation  2  must  be  made.  This  requires  the  replacements 

-*  c‘.\u:;i\(wu,+wm)  \ 

and  l  (10) 

C->  -  C-7  c*J  |(wv,+wv,).  j 

Analogous  substitutions  in  the  forward  directions  of  m  and  n  fol¬ 
low  automatically  in  the  integration  process. 
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The  velocity  replacements  in  Equations  5  may  be  illustrated  by 
the  example 


>  0,  C£V»,  >  0,  >  0,  H#.,  *  0, 

=0,?  + o . 

•»,»— 1 

One  finds  w  >  0,w  >  0,  and 

Vg-Uj*  (1+TPAM, 
U£„-*U£„(l-w), 

v&-*  v&a  +wb*a), 

and 


|  do 


i  (12) 


At  tiiis  point,  it  must  be  mentioned  that  attempts  were  explored 
to  lift  the  conrol  limits  prescribed  by  ki  and  ka  in  Equations  7b  and 
8b  in  an  effort  to  close  the  continuity  gap  completely.  However, 
since  the  bottom-friction  coefficient,  b,  is  rather  small  the  control 
limits,  ki  and  ka,  had  to  be  kept  small  to  achieve  best  results.  Com¬ 
putations  conducted  with  large  control  limits  k*  (excessive  bottom 
friction)  seemed  to  close  the  continuity  gap,  but  the  tidal  and  ve¬ 
locity  fields  in  the  open  oceans  assumed  unrealistically  small  values. 
Large  control  limits  ka  (insufficient  bottom  friction)  produced 
strong  instabilities  as  anticipated  from  the  analysis  in  the  section, 
“Stability  Analysis,’’  of  Part  L  To  safely  check  the  possible  insta¬ 
bility,  the  second  control  parameter  w  (Equations  5,  6,  and  12) 
was  defined  in  units  ofu  =  A4  andv  =  B\  in  contrast  to  u  =  1  and 
v-  1,  used  for  the  first  control  parameter  w. 

After  some  trial-and-error  computations,  the  following  control 
limits  were  chosen  for  die  ms-tide  model: 


ki  =  0.03,  ka  =  0.06. 


(13) 
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These  moderate  values  reflect  the  well-known  fact  that  the  magni¬ 
tude  of  bottom  friction  has  a  strong  effect  on  the  motions  consid¬ 
ered.  Indeed,  with  some  minor  improvements  of  the  tidal  field, 
significant  improvements  of  the  continuity  gap,  velocity  field,  and 
convergence  of  the  integration  were  achieved.  This  procedure  was 
applied  to  all  oceanic  cells  with  known  empirical  tide  data  (Equa¬ 
tion  3),  provided  these  cells  bordered  terrestrial  cells  or  contained 
small  islands  or  other  bottom  irregularities.  No  meaningful  reason  , 
was  seen  to  apply  the  same  bottom-friction  adjustment  procedure 
to  distant  offshore  oceanic  cells  with  available  deep-sea  tide  mea¬ 
surements. 

In  order  to  implement  the  second  step  (B)  of  the  hydrodynami- 
cal  interpolation  procedure,  the  following  velocity  replacements  in 
oceanic  mesh  cells  bordering  terrestrial  cells  were  defined: 


and 


Ul++\,  %U,U*X, 


\ 


>  (14) 


provided  in  Equation  3.  The  parameters  (uv)  are  mutually 
consistent  by  definition: 

Si  -  1  if  =  0,  otherwise  Si  =  0,  | 

2a  =  1  if  =  0,  otherwise  5»  =  0,  i 

Vi  =  1  if  *  0,  otherwise  Vi  =  0,  I  ^  ^ 

and 


Vi  =  1  if  =  0,  otherwise  =  0. 


The  remaining  continuity  gap  will  be  further  narrowed  when  die 
control  parameter  w  is  determined  to  be  in  agreement  with  Equa¬ 
tions  2, 4, 7, 8,  and  9  by 


% 


« 


V 

I' 
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{[ACtti-  "t  -  W{]/(  for  {Vo 

(16a) 

0  for  {  =  0 

with  the  third  control  limit 

\w\<k»,  (16b) 

where 

c = cu Sxittsu  -  u.t/£i]  +  vian«_i  - (i7) 

Obviously,  the  substitutions  (Equations  14)  specify  consistent 
in-  or  out-flows  across  the  mathematical  boundaries  of  oceanic 
coastal  cells,  as  illustrated  in  Figure  1,  without  explicitly  fixing  the 
physical  boundary  line.  Again,  no  complete  removal  of  the  con¬ 
tinuity  gap  was  possible.  The  most  satisfactory  results  for  the  Ms- 
tide  were  achieved  by  setting  the  third  control  limit  (Equation 
16b)  at 


ks  =  0.5:  (18) 

While  die  improvement  of  the  tidal  field  was  again  moderate,  the 
remaining  continuity  gaps  and  nearshore  velocity  distortions  as¬ 
sumed  uniformly  satisfactory  levels.  The  remaining  small  short¬ 
comings  of  die  model  can  easily  be  attributed  to  the  boundary 
inaccuracies  (c),  (d),  and  (e)  listed  above,  but  for  which  no 
simple  remedies  were  found. 

It  may  be  emphasized  that  the  rather  significant  change  in  die 
nearshore  velocity  field  permitted  by  the  in-  and  out-flow  specifica¬ 
tions  (Equations  14)  affected  the  tidal  field  only  in  a  minor  fashion. 
This  important  phenomenon  is  in  agreement  with  the  well-known 
fact  that  the  pressure  distribution  in  a  fluid  modem  is  very  insensi¬ 
tive  to  lane  but  local  velocity  variations.  For  instance,  it  is  perhaps 
the  most  important  postulate  in  Prandtl’s  boundary-layer  theory 
(see,  e.g.,  SchUchting,  1968  I),  and  it  is  the  basis  of  die  hydro- 
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static-pressure  assumption  invoked  here  and  in  the  section,  “The 
Continuous  Ocean-Tide  Equations  (COTEs),”  of  Part  I  for  the 
present  tidal  model. 

The  hydrodynamical  interpolation  technique  considerably  ac¬ 
celerated  the  convergence  of  the  integration  procedure  toward  the 
steady  state  amplitudes  and  phases.  In  fact,  the  computation  of  the 
new  M2-tide  model  (sections,  “Quality  of  the  Ocean-Tide  Model” 
and  “Conclusions”)  was  terminated  when  the  amplitudes  and 
phases  over  all  open  ocean  areas  differed  by  less  than  1  cm  and  1  °, 
respectively.  Obviously,  this  improved  convergence  feature  goes 
significantly  beyond  the  same  property  described  in  Part  I,  section, 
“Lateral-Boundary,  Initial  and  Final  Data,”  for  the  purely  mathe¬ 
matical  model. 

Quality  of  the  Ocean-Tide  Model 

Since  the  present  tide  model  incorporates  essentially  all  known 
empirical  data  by  hydrodynamical  interpolation  (preceding  sec¬ 
tion),  no  direct  comparison  of  observed  and  computed  data  is 
feasible.  Nevertheless,  a  comprehensive  appraisal  of  the  reality  of 
the  present  tide  model  is  possible  by  inspecting  the  quality  of 
hydrodynamical  interpolation;  i.e.,  by  evaluating  the  “smoothness” 
with  which  the  computed  tide  “accepts  or  rejects”  the  empirical 
tidal  data.  In  fact,  the  smoothness  characteristics  of  the  novel 
hydrodynamical  interpolation  technique  are  distinctly  different 
from  those  of  other  direct  interpolation  procedures  using  power  or 
trigonometric  polynomials.  In  the  latter  case,  smoothness  of  the 
interpolation  can  be  carried  up  to  any  desired  degree  by  simple 
design.  The  adjustment  of  hydrodynamical  parameters  (preceding 
section)  in  the  former  method  does  not  imply  any  smoothness  of 
the  interpolation,  unless  both  the  empirical  input  data  and  the 
hydrodynamical  tide  model  are  compatible  with  each  other. 

As  is  well  known,  local  tidal  distortions,  caused  by  an  isolated 
roughness  (seamount  or  small  island)  in  the  bottom  relief,  affect 
the  surrounding  ocean  tide  very  little.  The  major  level  of  ocean 
tides  is  shaped  by  continental  shorelines  and  large  (in  area  and/or 
length)  islands  and  ridges.  In  contrast  to  ordinary  polynomial  in- 
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terpolations,  an  important  feature  of  the  new  hydrodynamical 
interpolation  method  is  that  it  preserves  those  significant  properties 
of  ocean  tidal  currents  without  any  essential  alterations. 

Extensive  computer  experiments  were  conducted  to  test  the  im¬ 
portant  smoothness  characteristics  of  the  hydrodynamical  interpo¬ 
lation  procedure.  Faulty  input  data  were  deliberately  inserted  and 
quickly  recognized  as  rejected  by  the  computed  surrounding  tide. 
Indeed,  the  first  computations,  which  included  empirical  tidal  data, 
revested  immediately  several  input  errors  in  the  data.  Vice  versa, 
smoothly  accepted  empirical  tidal  data  were  randomly  deleted  to 
test  their  backlash  reaction  on  the  computed  tide.  As  anticipated, 
no  significant  modifications  were  detected.  Consequently,  the  hy¬ 
drodynamical  interpolation  technique  permits  a  check  of  the  reality 
of  both  the  tide  model  and  the  empirical  tidal  input  data.  If  an 
input  value  is  rejected  by  the  computed  tide,  then  one  or  the  other 
or  both  are  defective.  Fortunately,  only  very  few  discrepancies  be¬ 
tween  the  different  sources  of  observed  Ms-tide  data  (see  section, 
"Empirical  Tide  Data”)  have  been  discovered  that  way. 

The  new  discrete  tide-model  has  been  applied  to  compute  the 
global  Ms  ocean  tide.  A  complete  discussion  and  tabulation  of  all 
amplitudes  and  phases  is  presented  in  Schwiderski  (1979c  I).  In 
order  to  display  the  quality  of  die  tidal  model,  the  computed  am¬ 
plitudes  (in  cm)  and  phases  (in  degrees)  along  with  their  adjacent 
empirical  values  have  been  tabulated  in  ”30°  by  50°  map  form” 
for  four  typical  ocean  areas  (Tables  4-7).  All  empirically  sup¬ 
ported  input  data  along  continental  shores  and  at  island  stations 
are  underlined  in  the  tables.  All  nearshore  deep-sea  measurements 
deluded  in  the  model  are  labeled  by  subbrackets.  As  was  explained 
hi  die  preceding  section,  all  distant  offshore  deep-sea  measurements 
are  not  included  in  the  tide  model.  However,  their  approximate 
locations  are  marked  by  wavy  underlines,  and  their  corresponding 
observed  data  are  Haled  in  Table  3.  Land  points  are  left  blank. 

In  die  evaluation  of  the  tidal  accuracy,  one  must  remember  diet 
the  ocean  tide  at  any  fixed  location  is  determined  by  two  harmonic 
constants.  If  (fc,M  and  (1,1)  denote  the  respective  local  atnpli- 
todes  and  phases  of  the  "true”  and  "computed”  tides 
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b  =  (o  C®s(<rf  ~So),(  =  (  COS (vt  “  I), 


then  their  time-dependent  error  is 

?=io-:  =  fc08(<r/-I) 


with  die  standard  deviation 

rms(C)  =  V6V2I, 

where 


?  =  6.,-2^cos(*o-*)+f* 


and 


-  fe  sin  to  - I  sin  a 
““  Id  cos  Id  ~  I  cos  3  * 


Some  maximum  errors  are 


Car  *  Sr  ®  li  + 1  for  Id —  *  =  180°, 

&,  =  Sf  =lD-lfor«o-«  =  0°, 
tw  =  I*  =  2|  «n  V4  (Id  -  «)  for  I  =  Id, 


and 


{»  =  I*  =  I  for  I  —  Id  and  I*  —  •  =  60°. 


(19) 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

(27) 


Equation  27  expresses  die  important  (act  dial  a  60°  phase  error 
results  in  an  amplitude  error  equal  to  the  tidal  amplitude  and, 
hence,  renders  dm  computed  tidal  prediction  completely  useless.  Of 
course,  in  regions  of  sufficiently  small  amplitudes,  any  phase  error 
is  arceptaWs 

Tables  4A  and  4B  depict  the  tidal  amplitudes  and  phases,  re¬ 
spectively,  of  the  northwestern  Atlantic  Ocean  including  the  east¬ 
ern  Caribbean  Sat.  As  can  be  verified  by  earlier  tide  models,  this 
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entire  area  was  very  difficult  to  model,  because  its  rough  bottom 
topography  has  a  strong  effect  on  the  tidal  currents  that  sweep  over 
or  across  various  barriers  with  rapidly  changing  water  levels.  There 
is  the  broad  and  shallow  continental  shell  along  the  whole  North 
American  shoreline,  with  Cape  Hatteras,  Long  bland,  Cape  Cod, 

Nova  Scotia,  and  Newfoundland  all  protruding  into  the  ocean 
basin.  Furthermore,  there  are  the  Grand  Banks,  die  Bahama 
Banks,  and  the  long  and  narrow  Caribbean  Ridge.  Obviously,  all 
of  die  corresponding  local  tidal  features  could  not  be  realistically 
captured  by  the  tide  model  without  a  proper  representation  of  die 
bathymetry  (section,  MHydrodynamically  Defined  Ocean  Bathym¬ 
etry”)  and  without  the  hydrodynamkal  interpolation  (preceding 
section)  of  the  locally  collected  tidal  observations.  1 

Now,  if  one  scans  the  tidal  amplitudes  and  phases  (Tables  4A  \ 
and  4B)  from  the  north  to  the  south,  one  gathers  the  impression  f 

that  the  whole  computed  ocean  tide  is  completely  locked  into  the  f 

array  of  empirical  (underlined)  tidal  data  everywhere  along  die  j 

continental  coast  and  along  the  many  aligned  islands  separating  j 

die  Adantic  Ocean  from  the  Gulf  of  Mexico  and  die  Caribbean  i 

Sea.  It  is  particularly  impressive  to  see  the  observed  tide  data  at  ! 

the  offshore  islands  (Sable — SI,  Barbados — BB,  and  even  as  far  as  | 

Bermuda — BI)  and  at  die  included  nearshore  (subbrackets),  deep-  1 

sea  stations  all  realistically  well-accepted  by  die  computed  sur-  j 

rounding  tide.  Moreover,  one  finds  the  excluded  offshore  deep-sea 
measurements  (locations  marked  by  wavy  underlines)  in  the  At-  I 
lantic  and  Caribbean  Sea  fully  verified  by  the  independent  tide  1 
model. 

As  can  be  seen  in  the  special  listing  of  Table  3,  die  measured 
and  computed  amplitudes  and  phases  at  die  Adantic  stations  agree 
within  2  cm  and  6°,  respectively.  The  remaining  discrepancy  is 
probably  within  die  experimental  error  due  to  short  observation 
times  and  die  use  of  die  distant  reference  station  Bermuda  (Zettior 
tt  at.,  1975),  which  exhibits  even  larger  gaps  between  the  various 
tidal  observations  listed  in  Table  2. 

Attention  may  be  drawn  to  die  existence  of  considerable  slopes 
between  the  empirical  boundary  data  and  the  computed  ocean-tide 
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values  in  die  high-amplitude  ranges  from  Nova  Scotia  to  Cape  Cod 
and  from  Cape  Hatteras  to  Florida’s  coast.  Yet,  these  rapid  tidal 
variations  can  be  considered  as  realistic  because  throughout  the 
same  sections  the  empirical  data,  among  themselves,  display  ex¬ 
actly  the  same  roughness.  This  only  substantiates  clearly  the  fun¬ 
damental  difference  between  polynomial  and  hydrodynamical  in¬ 
terpolation  techniques  pointed  out  above. 

In  the  complete  report  (Schwiderski,  1979c  I),  the  same  tidal 
roughness  will  be  recognized  in  several  similar  coastal  places  around 
the  world.  From  this  typical  phenomenon,  one  can  draw  the  fortunate 
conclusion  that,  while  some  empirical  data  may  be  lacking  high 
accuracy  (see  Table  1  and  the  British  Admiralty  Tide  Tables, 
1977  I),  the  computed  adjacent  ocean  tide  may  retain  its  high 
quality. 

In  order  to  gain  a  deeper  insight  into  the  detailed  tidal  phenom¬ 
ena  from  the  enclosed  table  charts  (e.g..  Tables  4A  and  4B),  it  is 
helpful  to  recall  the  physical  meaning  of  the  tabulated  tidal  con¬ 
stants.  The  local  tidal  amplitude,  is  defined  as  half  the  tidal 
“range,”  which  measures  the  total  variation  of  the  water  level  from 
high  to  low.  Lines  Of  constant  amplitudes  are  called  “corange 
lines.”  The  local  phase,  I,  specifies  the  tidal  cresting  time  (in  de¬ 
grees)  after  the  moon’s  (or  sun’s)  passage  over  the  Greenwich 
meridian.  For  the  present  Ms-tide  one  has  the  following  time 
conversions: 

360°  =  12.421  hr  (period), 

30*=  1.035  hr,  (28) 

1°  =  2.070  min. 

Lines  of  constant  phases  (simultaneous  cresting  times)  are  called 
“cotidal  lines.”  In  particular,  at  the  0°  *  360*  cotidal  lines,  which 
are  conspicuously  viable  in  the  phase  charts  (Tables  4B  to  7B), 
the  tide  crests  simultaneously  with  the  moon’s  passage  over  the 
Greenwich  meridian.  The  tidal  crest  advances  with  time  normal  to 
die  cotidal  lines  toward  larger  phases.  A  point  of  aero  amplitude 
(f  18  0)  around  which  the  tidal  crest  rotates  from  0*  to  360°  is 
called  an  “amphidromic  point”;  it  is  marked  in  the  tables  by  a 
circled  star  (g) . 
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Key  to  M2-tide  Tables  4-7. 

=  Longitude  east  (°). 

=  Colatitude  (°). 

=  Amphidromk  point. 

=  Subbars  mark  empirical  input  data. 

=  Subbrackets  mark  input  nearshore  deep-sea  measurements. 

=  Wavy  underlines  mark  offshore  deep-sea  tide  gauge  stations  with 
excluded  measurements  listed  in  Table  3a,  b. 
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Ocean  Tides,  Part  II:  A  Hydrodynamics!  Interpolation  Model 
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In  the  area  of  Tables  4A  and  4B,  a  major  amphidromic  point  is 

visible  in  the  Caribbean  Sea  southeast  of  the  island  of  Puerto  Rico 

•* 

(PRI)  near  the  marked  deep-sea  gauge  station.  The  loosely  con¬ 
nected  Caribbean  and  Atlantic  tides  rotate  counterclockwise 
around  this  point  with  the  0°  =  360°  cotidal  line  running  north¬ 
eastward.  As  a  result  of  this  rotation,  the  whole  Caribbean  Sea 
appears  to  be  trapped  and  unable  to  develop  any  significant  Ms- 
tide.  In  agreement  with  the  observations,  the  Ms  tidal  crest  sweeps 
across  the  Caribbean  Sea  essentially  from  north  to  south  with  very 
little  variation  in  water  level. 

If  one  follows  the  tidal  crest  around  the  amphidromic  point 
from  the  Atlantic  Ocean  to  the  Caribbean  Sea  and  back  to  the 
Atlantic,  one  recognizes  a  major  tidal  distortion  caused  by  ocean 
ridges,  which  has  long  been  discovered  by  practical  tidalists  (see, 
e.g.,  Harris,  1904 1;  Bogdanov,  1961 1;  Defant,  1961 1;  and  Luther 
and  Wunsch,  1974  I).  As  the  tide  crosses  the  ridge  between  the 
islands,  it  suffers  a  distinct  amplitude  jump  and  a  significant  phase 
shift  For  example,  north  of  Puerto  Rico  (PRI)  and  Hispaniola 
and  in  die  southeast  around  Barbados  (BB),  the  computed  and 
empirical  Atlantic  tide  data  display  a  higher  water  level  and  an 
earlier  or,  respectively,  a  delayed  cresting  time  than  the  adjacent 
tide  data  on  the  Caribbean  side.  In  particular,  in  full  agreement 
with  the  observations,  the  tidal  retardation  time  can  easily  exceed 
30°  (~1  hour).  The  distortion  seems  to  depend  on  the  angle  with 
which  die  tidal  crest  spills  over  the  ridge.  Maximum  distortion  ap¬ 
pears  to  be  associated  with  a  normal  crossing.  It  may  be  pointed 
out  that  the  realistic  resolution  of  tidal  distortions  by  ocean  ridges 
(see  below)  constitutes  probably  the  most  significant  improve¬ 
ment  of  die  new  model  over  all  earlier  hydrodynamical  models. 

The  Atlantic  portion  of  the  Caribbean-Atiantic  amphidromic 
rotation  is  opposed  by  a  southward  advancing  tide  Bom  about 
Newfoundland  in  die  north  and  by  an  eastward  progressing  tide 
from  about  Cape  Cod  to  Cape  Hatteras  in  the  west  As  a  result  of 
this  interaction  of  three  opposing  tidal  waves,  the  middle  latitudes 
(around  n  —  60°)  of  the  Atlantic  display  very  small  variations  in 
tidal  amplitudes  and  phases.  In  the  high-amplitude  sections  be- 
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tween  Nova  Scotia  and  Cape  Cod  and  between  Cape  Hatteras  and 
Florida’s  coast,  the  Caribbean-Atlantic  rotation  wave  seems  to  be 
less  affected  by  the  opposing  tidal  waves  and  progresses  frontally 
against  the  corresponding  shallow  coastal  corners. 

Since  the  tide-generating  Ms-potential  is  a  single  progressing 
wave  from  east  to  west,  the  ocean  responds  with  amphidromic  tidal 
waves  that  cannot  reverse  their  directions.  Thus,  at  shore  points 
tidal  waves  are  either  incoming  or  outgoing  without  reversals.  In 
the  first  case  tidal  crests  always  move  from  sea  to  shore.  In  the  sec¬ 
ond  case  tides  always  swell  to  their  crests  at  the  shore  first  and  then 
move  out  to  sea.  The  incoming  tide  between  Nova  Scotia  and  Cape 
Cod  seems  to  produce  high  and  rough  waters.  The  outgoing  tide 
between  Cape  Cod  and  Cape  Hatteras  is  distinctly  lower. 

Although  the  computed  tide  in  the  Gulf  of  St.  Lawrence  displays 
the  well-known  amphidromic  point  (Defant,  1961  I),  the  grid  sys¬ 
tem  is  much  too  crude  to  attach  a  high  accuracy  to  the  tidal  con¬ 
stants  in  this  border  sea.  For  the  same  reason,  the  tidal  data  listed 
between  Florida,  Cuba,  and  the  Bahamas  are  naturally  less  accu¬ 
rate  than  those  in  the  open  oceans. 

Tables  5A  and  SB  illustrate  the  smoothness  with  which  die  com¬ 
puted  tide  of  the  northeastern  Pacific  Ocean  attaches  itself  to  the 
empirical  tide  data  along  the  North  American  west  coast  The  tidal 
constants  observed  at  the  islands  of  Guadalupe  (GI)  and  Farallon 
(FI),  at  the  Cobb  Seamount  (CS),  and  at  the  included  nearshore 
deep-sea  stations  fit  realistically  well  into  the  computed  surround¬ 
ing  tide.  The  amplitudes  and  phases  of  the  excluded  offshore  deep- 
sea  measurements  in  die  Pacific  agree  within  2  cm  and  6°,  respec¬ 
tively,  with  the  computed  data  (Table  3),  which  is  just  the  same 
accuracy  as  in  die  Atlantic. 

Perhaps  die  most  prominent  feature  of  this  area  is  the  amphi¬ 
dromic  point  ® ,  around  which  the  Ms-tide  rotates  counterclock¬ 
wise.  This  amphidromic  system  was  predicted  by  Munk  el  al. 
(1970)  and  Irish  et  al.  (1971)  in  almost  identical  geographical 
position.  Earlier  hydrodynamical  tide  models  failed  to  resolve  this 
system  on  proper  location,  although  several  models  matched  the 
empirical  data  along  die  coast  quite  well.  Since  die  northeastern 
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Pacific  falls  short  in  major  bottom  and  coastal  irregularities  when 
compared  to  the  northwestern  Atlantic,  the  indicated  rapid  loss  of 
quality  in  westerly  direction  seemed  disappointing.  Yet,  as  will  be 
demonstrated  below,  this  shortcoming  could  have  been  concluded 
from  the  obvious  failure  of  those  models  to  reasonably  reproduce 
the  tide  over  most  of  the  north  and  central  Pacific  Ocean. 

As  was  mentioned  before,  the  author’s  preliminary  tide  model 
(Schwiderski,  1976  I)  used  a  bathymetry  that  failed  to  represent 
the  hydrodynamical  barrier  effects  of  the  Marianas,  Nampo,  Kuril, 
Aleutian,  and  Hawaiian  ridges,  as  well  as  of  other  seamount  chains. 
Consequently,  the  Mi -tide  of  almost  the  whole  central,  western, 
and  northern  Pacific  area  was  modeled  as  a  single  huge  amphidro- 
mic  system,  as  pictured  by  the  similar  maps  of  other  numerical 
tidalists  such  as  Zahel  (1971  I)  and  Estes  (1975  I,  1977  I).  The 
clockwise-rotating  Pacific  tide  was  free  to  sweep  undisturbed  into 
the  Philippine,  Okhotsk,  and  Bering  seas.  By  the  time  the  computed 
tidal  crest  reached  the  Aleutian  Islands,  it  was  just  about  180°  out 
of  phase.  When  the  original  bathymetry  was  replaced  by  hydrody- 
namically  defined  depth  data  (section,  “Hydrodynamically  Defined 
Ocean  Bathymetry”),  the  entire  Pacific  Ocean  resembled  a  whirl¬ 
pool  after  some  continued  computations  over  several  quarter  pe¬ 
riods.  The  amphidromic  system  weakened,  and  its  center  slipped 
slowly  southward,  but  drastically  improved  phases  appeared  grad¬ 
ually  along  the  Aleutian  Ridge,  confirming  the  anticipated  effect  of 
ocean  ridges. 

The  complete  turnaround  of  the  Pacific  Ms  tide  near  the  Aleu¬ 
tian  Islands  was  speeded  up  when  the  empirical  tidal  constants  were 
introduced  into  the  model.  In  fact,  a  repeat  of  the  same  computa¬ 
tions  settled  the  Pacific  Ocean  tide  into  its  final  position  in  a  rather 
dramatic  fashion.  Striking  improvements  woe  registered  over  the 
whole  Pacific  and,  of  course,  also  over  the  Atlantic  and  Indian 
oceans. 

As  is  depicted  in  Tables  6A  and  6B  for  the  north-central  Pacific, 
the  amphidromic  system  is  replaced  by  a  low-amplitude  tide.  It 
appears  to  be  locked  in  between  the  Aleutian  and  Hawaiian  ridges 
in  the  north  and  south  and  also  between  the  Emperor  Seamount 
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chain  in  the  west  and  the  high-amplitude  tide  in  the  east,  which 
progresses  in  a  westerly  direction  from  the  west  coast  of  North 
America  (Tables  SA  and  SB).  The  amplitude  topography  of  this 
area  resembles  the  low-amplitude  tide  in  the  Caribbean  Sea  (Table 
4A).  When  the  westward-advancing  tidal  wave  enters  the  region 
between  the  Aleutian  and  Hawaiian  ridges,  it  suffers  a  remarkable, 
almost  symmetric  retardation  at  both  ridges.  In  fact,  as  the  visible 
(0°  =  360°)  cotidal  line  in  Table  6B  reveals,  the  crest  front  of  the 
tidal  wave  assumes  the  shape  of  an  almost  symmetric  wedge.  If  one 
traces  the  0°  phase  line  westward  beginning  at  both  ridges,  one  can 
infer  a  definite  idea  about  the  realistic  reproduction  of  the  tide  in 
this  region.  At  both  ends,  the  0°  phase  is  in  full  agreement  with  die 
empirical  data.  As  the  observed  phases  grow  westward  along  both 
ridges,  so  grow  proportionally  the  distances  of  the  0°  phase  line 
from  the  ridges. 

The  new  computed  Ms-tide  model  no  longer  indicates  any  symp¬ 
toms  of  the  original  phase  problems  at  the  Aleutian  and  Hawaiian 
ridges.  The  computed  amplitudes  and  phases  approach  the  empiri¬ 
cal  tidal  constants  from  both  sides  of  the  ridges  as  smoothly  as 
could  be  desired.  As  the  tidal  wave  spills  over  both  ridges  in  north¬ 
westward  or  southwestward  directions,  respectively,  it  suffers  a 
tidal  distortion  similar  to  that  found  before  at  the  Caribbean  Ridge. 
Amplitude  jumps  and  major  phase  shifts  are  again  in  complete 
agreement  with  observations  (see  the  remarks  of  Luther  and 
Wunsch,  1974 1).  It  is  particularly  gratifying  to  find  the  phase  shift 
well  developed  along  the  whole  length  of  the  Hawaiian  Ridge  from 
the  island  of  Hawaii  to  Midway,  even  though  only  few  stations  of 
data  were  used  at  both  ends.  Also,  it  may  be  noticed  that  the  ob¬ 
served  tidal  constants  at  die  distant  and  isolated  island  stations  of 
Pribilof  (PF),  Midway  (MW),  and  Johnston  (JI)  are  all  realisti¬ 
cally  well  integrated  by  the  surrounding  computed  tide. 

Ironically,  the  old  and  new  Ma-tide  maps  constructed  by  Bog¬ 
danov  (1961 1)  and  Luther  and  Wunsch  (1974  I)  by  pure  intui¬ 
tion  and  simple  rules  of  thumb  from  empirical  data  came  closest  to 
die  present  charts.  Indeed,  their  maps  display  no  amphidromic 
system  in  the  north-central  Pacific.  As  is  verified  in  Schwiderski 
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(1979c  I),  the  computed  amphidromic  points  between  the  Cook 
and  Society  islands  and  near  the  southern  edge  of  the  Solomon 
Islands  are  both  in  almost  identical  positions  with  those  charted  by 
the  same  authors.  Nevertheless,  their  detailed  distribution  of  ampli¬ 
tudes  and  phases  is  still  significantly  different  from  the  present  one. 

Perhaps  the  most  spectacular  display  of  the  high  quality  of  both 
the  computed  and  the  observed  tidal  data  is  brought  out  by  Tables 
7A  and  7B  depicting  the  high-amplitude  tide  of  the  central  Pacific. 
Indeed,  unlike  any  other  open  ocean  area,  the  tabulated  region  is 
dotted  with  numerous  tide-gauge  stations  at  island  groups  and  at 
scattered  isolated  islands.  In  addition  to  the  fully  listed  island 
chains,  there  are  the  isolated  islands:  Johnston  (JI),  Wake  (WI), 
Kudaie  (KI),  Ocean  (OI),  Funafuti  (FI),  Wallis  (WI),  Niue 
(NI),  and  Norfolk  (NF).  The  corresponding  observed  tidal  con¬ 
stants  listed  in  nongeographical  arrangement  appeared  incoherent 
and,  hence,  uncorrelated,  giving  rise  to  doubt  their  true  value.  Yet, 
the  computed  tidal  wave  sweeps  across  the  whole  area  in  a  south¬ 
westerly  direction  with  little  variation  of  its  high  amplitude.  As  the 
wave  crest  passes  through  the  many  checkpoints  with  correct  height 
and  in  right  time,  it  integrates  and  correlates  without  a  single  ex¬ 
ception  all  the  empirical  data  into  one  coherent  unity. 

Conclusions 

The  quality  evaluation  of  die  constructed  Mi  ocean  tide  model 
described  in  Parts  I  and  II  of  this  paper  leads  to  the  conclusion  that 
it  is  now  possible  to  compute  detailed  and  accurate  global  ocean 
tides  which  fulfill  die  application  requirements  of  contemporary 
researchers.  In  fact,  it  is  estimated  that  the  computed  Mi-tide  charts 
permit  an  Mi-tide  prediction  anywhere  in  the  open  oceans  with 
an  accuracy  of  better  than  5  cm.  This  accuracy  leaves  ample  room 
for  superposable  errors  due  to  the  additional  smaller  tidal  constitu¬ 
ents  listed  in  Table  1  of  Part  I,  which  are  presently  under  construc¬ 
tion  with  equivalent  relative  accuracy.  When  all  those  partial  tides 
become  available,  the  total  tide-prediction  error  is  expected  to  fall 
well  below  the  10-cm  limit  needed  in  many  applications. 

Naturally,  the  achieved  high  accuracy  of  the  Mt- tide  in  the  open 
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oceans  drops  somewhat  near  continental  or  island  stations  where 
empirical  data  are  missing  or  are  less  accurate  themselves  (see  the 
introduction  to  the  British  Admiralty  Tide  Tables,  1977  I).  This 
is  particularly  true  near  Antarctica  and  in  the  Arctic  Ocean,  where 
reliable  measurements  of  ocean  tides  and  depths  are  sparse.  Also, 
less  accurate  predictions  must  be  anticipated  in  small  border  seas, 
bays,  estuaries,  and  channels  where  the  1°  by  1°  grid  system  pre¬ 
cludes  a  sufficient  resolution.  To  improve  the  present  tide  model  in 
those  areas,  significantly  improved  observations  will  be  needed 
along  with  a  locally  refined  network  and  corresponding  bathy¬ 
metric  data. 
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